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Abstract 

This work concerns single-trace correlations of Euclidean multi-matrix models. In 
the large- N limit we show that Schwinger-Dyson equations (SDE) imply loop equa- 
tions (LE) and non-anomalous Ward identities (WI). LE are associated to generic 
infinitesimal changes of matrix variables (vector fields). WI correspond to vector 
fields preserving measure and action. The former are analogous to Makeenko-Migdal 
equations and the latter to Slavnov- Taylor identities. LE correspond to leading 
large- ./V SDE. WI correspond to 1/N 2 suppressed SDE. But they become leading 
equations since LE for non-anomalous vector fields are vacuous. We show that sym- 
metries at N = oo persist at finite N , preventing mixing with multi-trace correla- 
tions. For 1 matrix, there are no non-anomalous infinitesimal symmetries. For 2 or 
more matrices, measure preserving vector fields form an infinite dimensional graded 
Lie algebra, and non-anomalous action preserving ones a subalgebra. For Gaussian, 
Chern-Simons and Yang-Mills models we identify up to cubic non-anomalous vec- 
tor fields, though they can be arbitrarily non-linear. WI are homogeneous linear 
equations. We use them with the LE to determine some correlations of these mod- 
els. WI alleviate the underdeterminacy of LE. Non-anomalous symmetries give a 
naturalness-type explanation for why several linear combinations of correlations in 
these models vanish. 
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1 Introduction 



Hermitian multi-matrix models are quantum systems where the dynamical variables are a set 
of N x N hermitian matrices. Observables must be basis independent, i.e. invariant under the 
global adjoint action of U(N) on the matrices. Expectation values of observables are determined 
by an average over all matrix elements with respect to a Boltzmann weight specified by an 
action. Matrix models simplify in 't Hooft's large- N limit, since fluctuations in U(N) -invariant 
observables are small in this limit. 

Multi-matrix models are simplified models for the dynamics of gauge fields in Yang-Mills 
theory. It is a fundamental and challenging problem to determine the free energy and correlation 
functions of multi-matrix models, and elucidate the mathematical framework needed to study 
them. Multi-matrix models are much harder to understand than single-matrix models, but 
also have a much richer structure. Large- N matrix models and more generally large- N gauge 
theories have been studied ever since their relevance as an approximation to the theory of strong 
interactions was pointed out by 't Hooft in the mid 1970s [U [21 [3]. Important progress in 
obtaining the loop equations of Yang- Mills theory and study of the large- N limit was made in 
the late 1970s and early 1980s by Migdal and Makeenko [3J, Cvitanovic[5], Yaffe[6], Jevicki and 
Sakita[7J and others [8]. The subject was applied to random surface theory, 2d string theory and 
the matrix approach to M-theory in the 1990s. Meanwhile, there has been a steady stream of 
developments in matrix models of which we cite a few examples. These include their connections 
to non-commutative probability theory [THl [TT1 [T2J [13l [TJ] , the study of multi-matrix symmetry 
algebras and their connections to spin chains [151 116j. exact solutions [17] and their relation to 
CFT[18j and algebraic geometry and detailed studies of the loop equations [191 [20] . Much of 
the existing literature deals with 1-matrix models or exact solutions for specific observables of 
carefully chosen multi-matrix models. We hope to complement this by developing a framework 
and methods that apply to general multi-matrix models. 

Throughout physics, we exploit symmetries to simplify dynamical equations by reducing the 
number of unknowns. The quantum dynamical equations of a large- N multi-matrix model are 
the loop equations^ for single-trace correlations. These correlations are analogs of gluon and 
ghost correlation functions of Yang-Mills theory. Here, we develop a general framework to find 
non-anomalous infinitesimal symmetries of multi-matrix models (i.e. those that preserve both 
action and measure). These symmetries are used to infer Ward identities, which supplement the 
loop equations to determine correlations. These non- anomalous symmetries and Ward identities 
can be regarded as finite dimensional analogues of BRST invariance and Slavnov- Taylor identities 
of Yang-Mills theory. The ideas are illustrated with examples from 2 and 3-matrix models. 

To motivate this work, we explain how we came to think along these lines. We were trying 
to solve the loop equations (LE) to determine large- N single trace correlations of some specific 
multi- matrix models [20]. Single-trace correlations are the basic objects of interest, since multi- 
trace correlations factorize into products of single-trace ones in the large- N limit. The LE state 
the invariance of the partition function under infinitesimal but non-linear changes of integration 
variables, in the large- N limit. They relate a change in action to a change in measure. Such 
infinitesimal changes of integration variables can be regarded as vector fields. A priori, there 

lr rhe name loop equations is used because these equations are analogous to the Makeenko-Migdal equations 
of Yang-Mills theory, which were formulated for the Wilson loops. The word loop has nothing to do with loops 
in Feynman diagrams. Another name for these equations is factorized Schwinger-Dyson equations. 
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is one LE for each such vector field. In many non-Gaussian cases, we found the LE were 
underdetermined. Moreover, this underdeterminacy seemed related to the fact that for several 
changes of variable, the LE were vacuous. In other words, for some vector fields, both the change 
in action and change in measure simultaneously vanished in the large- N limit. We were looking 
for additional equations to supplement the LE and alleviate their underdeterminacy. 

Now, the LE can be regarded as the large- N limit of the finite- N Schwinger-Dyson equations 
(SDE). The SDE are conditions for the invariance of matrix integrals for multi-trace correlations 
under infinitesimal non- linear changes of integration variables. They relate a change in action to 
a sum of a change in measure and change in observable. While the first two are usually of order 
iV , the latteid is usually of order N~ 2 . So naively, in the large- N limit, the latter drops out and 
using factorization, we get back the LeH. However, in the special case where the vector field is a 
symmetry of both action and measure in the large- N limit, the naive large- N limit of the SDE 
is vacuous and one must go to the next order in 1/N 2 to get the leading condition. For generic 
vector fields, this 0(1/ N 2 ) SDE would not be an equation for the single-trace correlations 
alone, since it involves 0(1/ N 2 ) corrections to the factorized result for multi-trace correlations. 
However, remarkably, we show that if a vector field is a symmetry of the action in the large- N 
limit, then it is also a symmetr>Q at each order in 1/N 2 . The same also holds for symmetries 
of measure. The simple reason is that there are more independent variables as ./V grows, and 
so more conditions on a vector field to be a symmetry as N increases. Non-anomalous vector 
fields define simultaneous symmetries of action and measure in the large- N limit. Thus, for 
non-anomalous vector fields, the change in action and measure terms drop out of the 0(1/N 2 ) 
SDE, which then becomes a condition for invariance of the expectation value of the observable in 
the large- N limit, schematically L v Gi x ...i n = 0. This latter condition only involves single-trace 
correlations Gi x ...i n , and is what we call a non-anomalous Ward identity (WI). WI are associated 
to vector fields v that leave both action and measure invariant in the large- N limit. L v is the 
Lie derivative along v. It is precisely for such vector fields that the LE are vacuous. We show 
that such vector fields form an infinite dimensional graded Lie algebra. Thus, as N — > oo, the 
SDE imply not just the LE, which are associated to generic vector fields, but also WI, which are 
associated to non- anomalous vector fields. The latter are easily overlooked in a naive passage 
to the large- iV limit. Naively, the WI appear to be 'universal', i.e. independent of the matrix 
model action, but this is not really true. Whether or not a WI holds is determined by whether 
the corresponding vector field is a simultaneous symmetry of both action and measure. 

In retrospect, these non-anomalous symmetries and WI are not unexpected. In Yang-Mills 
theory, non-anomalous symmetries include Poincare and BRST invariance. WI for the latter 
are Slavnov- Taylor identities. Our non-anomalous WI share with the Slavnov- Taylor identities 
the structural similarity of being homogeneous linear equations for correlations. Just like our 
WI, the Slavnov- Taylor identities seem independent of the gauge-fixed Yang-Mills action, until 
one realizes they hold only because the action and measure are BRST invariant. In Yang-Mills 
theory, while Poincare transformations act linearly on the fields, BRST transformations are 
quadratically non-linear. For specific matrix models, we find non- anomalous symmetries that 
are linear, quadratic (n = 2) and cubic (n = 3) ; there is no limit to the possible non-linearity 
of such symmetries. Moreover, for n > 1, some rank-n + 1 non-anomalous symmetries can be 

2 There is no 0(1/ N) contribution in a matrix model, i.e. in the absence of quarks or N -vectors. 
3 This is why LE are also called factorized Schwinger-Dyson equations. The name Virasoro constraints is also 
used, especially in applications to string theory. 

4 There are potentially more symmetries at finite- N than at large- N , so the converse is almost certainly false. 
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obtained via the Lie brackets of rank- re symmetries. This is reminiscent of how Poisson brackets 
of conserved charges (if non- vanishing) , give higher conserved charges in integrable models. 

We find a significant difference between 1 -matrix models and multi-matrix models. The 
measure for a single matrix in the large- N limit admits only one continuous symmetry, i.e. 
translations of the matrix. Translations, however, are not a symmetry of any non-trivial 1- 
matrix action. Thus, non-trivial 1 -matrix models have no non-anomalous WI. Interestingly, we 
find that the measures for multi-matrix models allow for large classes of symmetries, some of 
which may also be symmetries of a given action. 

In practice, once a non-anomalous symmetry of a model is known, it is easier to first solve 
the resulting WI and then consider the LE. The WI, being homogeneous linear equations, force 
several correlations or linear combinations thereof to vanish. This simplifies analysis of the 
LE, which are mildly non-linear. However, we caution that some WI may contain the same 
information as contained in the LE, while others may provide new conditions. 

We emphasize that the techniques and results of this paper are exact. They do not involve any 
approximation beyond the passage to N = oo. Our methods apply to single-trace correlations 
of generic hermitian multi-matrix models with polynomial actions. They are not special to any 
subclass of actions or correlations. Of course, the non-anomalous symmetries and WI will depend 
on which model we consider. One lesson we learned is that though it is a bit laborious, it is 
possible to solve the LE and WI of large- ./V multi- matrix models to determine exact correlations, 
starting from the lowest rank ones. 

Finally, our derivation of the SDE, LE and WI makes use of the matrix integral representation 
for correlations. In cases where these integrals converge, we expect the equations to be rigorously 
valid. When the matrix integrals diverge, the SDE, LE and WI are only formal statements and 
their consistency is not guaranteed by our work. Indeed, we seem to find an example where formal 
use of these equations for a model whose matrix integrals diverge, leads to inconsistencies. We 
are yet to understand the deeper significance of this. 

We can give another interpretation of our results. Suppose one were to calculate single trace 
correlations of a large- N multi-matrix model. Then in many cases one would find there are 
several linear combinations of correlations that vanish. One might look for a naturalness-type 
explanation for this i.e., a non-anomalous symmetry that forces those linear combinations to 
vanish. In some cases, there is a discrete symmetry (such as A — > —A for correlations of odd 
order for an even action) that does the job. The results of this paper may be regarded as the 
discovery of several new continuous non-anomalous symmetries of multi-matrix models. For 
example, 5A\ = a[A\, [^^l]] an d 5A2 = a[[Ai, A2], A2] is a non-anomalous symmetry of the 
Gaussian+YM 2-matrix model for all real a. Such symmetries lead to WI, which ensure that 
the quantities in question vanish. 

Organization and summary of results: In section [2] we determine the SDE of hermitian 
multi-matrix models. We show that in the large- N limit, they lead to LE supplemented by WI. 
It is asserted that the WI are to be imposed for every vector field that is a simultaneous symmetry 
of action and measure in the large- N limit. The proof of validity of the WI is completed in 
sections l3.3l and r4.11 In section l2.1l we explain why the WI trivialize for a 1 -matrix model. Section 
12.21 exhibits that multi-matrix LE are often underdetermined and this motivates the need for 
additional equations to determine correlations. WI potentially alleviate the under determinacy 
of LE. In section [3] we characterize measure preserving vector fields of multi-matrix models in 
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the large- N limit. We show that they form an infinite dimensional Lie algebra (section I3.4[) . 
Measure preserving transformations of 2- and 3 -matrix models are given in section [331 We work 
out the linear and quadratic non-anomalous symmetries of Gaussian, Chern-Simons, Yang-Mills 
and Gaussian-I- Yang- Mills multi-matrix models in sections 14.31 and 14.41 and also construct some 
cubic symmetries via Lie brackets of quadratic ones. In section we explicitly give the LE 
and non-anomalous WI for the Gaussian, Gaussian-I- YM and Chern-Simons models. We show 
that several correlations vanish, determine some non-vanishing correlations, and also obtain 
non-trivial relations among other non- vanishing correlations. In section [5~H we show that formal 
use of LE and WI for a model whose matrix integrals do not converge potentially leads to 
inconsistencies. Some outstanding questions are collected in [6j In appendix [A] we give an 
alternate derivation of the SDE that preserves hermiticity of matrices. In appendix IA.ll we 
consider some other possible changes of variables in an unsuccessful search for equations satisfied 
by the N = oo single trace correlations, over and above the LE and WI. In appendix[C]we argue 
that the WI by themselves (without use of LE) cannot determine all correlations of a non-trivial 
model. In appendix [D] we quote a useful formula for the number of cyclically symmetric tensors 
of rank-n in a A-matrix model. 



2 Schwinger-Dyson equations and Ward identities 

We consider a bosonicH Euclidean matrix model with A random hermitian matrices A{,i = 
1,2, • • • , A. The action tr S(A) = tr J2\j\< m S J Aj is taken to be a polynomial. Due to the 
trace, only the cyclic projections of the coupling tensors S 1 contribute to the action. Multi- 
indices are denoted by capital letters, for example, / = i\%2 ■■■in- Repeated lower and upper 
indices as summed and |/| denotes the length of the multi- index. 

Observables are functions of Ai that are invariant under the global adjoint action Ai — > 
UAiU^ of U € U(N). An important class of such functions are the trace invariants = 
-k tr Aj . The partition function and multi-trace correlations are defined as 

Z = J J] dA je ~ N tr ^ and • • • = I J UjdAje^ tr *W* Kl (1) 

(&Ki ■ * * &K n ) i s symmetric under interchange of any pair from K\, ■ ■ ■ , K n . It is cyclically 
symmetric in each Ki separately. (^k 1 ■ • • ®K n ) may be expanded in inverse powers of iV 2 

(* Kl ■ ■ ■ **J = G£U,.. ; *„ + jr 2 G% K2 .., Kn + j^G%i )K2 ^ Kn + ■■-. (2) 

The coefficient of N~ 2h can be regarded as a sum of Feynman diagrams that can be drawn on a 
Riemann surface with h handles and n disks cut ou®. The perimeter of each disk is associated 
to one of the inserted K^s. In particular, for h = 0, these are planar diagrams. Each of the 
Gg^)„.. K is symmetric in the multi-indices K\, • • • , K n . Factorization of multi-trace correlations 
in the large- A limit [22] means that G(°) can be written as a product of single trace correlations 

lim ($ Kl • • • $ Kn ) = G%\. K2 ..... Kn = G Kl --- G Kn , where G K = lim &*> K ). (3) 

5 It is possible to extend these methods to models with ghion and ghost fields, such as gauge fixed Yang-Mills 
theory. In this case the A, would include hermitian matrices as well as matrices with grassmann entries. 

6 A pictorial representation of ([2]) for a four-point correlation would resemble Fig. 1.8 on page 31 of Ref. [21] 
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The single-trace gluon correlations Gk are cyclically symmetric in K and satisfy the hermiticity 
condition G* K = G K provided S 1 also satisfy this property. Here K is the word K with order 
of indices reversed. Gk will also be referred to as moments, they are the moments of a non- 
commutative probability distribution [13] when A > 1. The rank of Gk is defined as \K\. 

To determine correlations, we derive Schwinger-Dyson equations (SDE), conditions for invari- 
ance of matrix integrals for (<&k x • • • *&K n ) under infinitesimal non- linear changes of variables^ 

[Ai\l -» [A'^l = [Ai\l + v{ [Ai\l, v\ infinitesimal real parameters for |/| > 0. (4) 

These include the BRST-type of transformations used to derive the Slavnov- Taylor identities of 
gauge-fixed Yang-Mills theory. For example, in Lorentz gauge the BRST transformations are 
infinitesimal quadratic transformations (A is an infinitesimal anti-commuting parameter), 

An^An + XdpC + XlA^c], c^c + A[c,c]+, c-*c + ASM M . (5) 

To calculate the effect of (jSJ) on the integral JTJ defining ($>Ki ' ' ' ^K n ) we neecHthe infinitesimal 
change in action, measure and the inserted observable $k 

e -NtrS'Aj e- NtisJAj (l-N 2 viS J ^j lIj2 ) + 0(v 2 ), 

det(|f||) = l + iVV^^^ + O^ 2 ) 

$K -» *K + 8g M vl* LIM . (6) 

The conditions for invariance of (&Ki ■ • • ®K n ) to linear order in vj are the finite N SDrfl 

v i S ^(<S> JlIj2 <S> Kl ---<S>K n ) = v!6 I / iI *(<f> h <!> l2 $ Kl ---$ Kn ) 

1 n 

+^2 E • ■ ■ *k p -Al p im p *k v+1 ■ ■ ■ *Jf„>.(7) 

p=l 

There is a priori one such SDE for each vector field v\ and each n = 0, 1, 2, where n is 
the number of insertions. The LHS is the expectation value of the change in action (along with 
$>K insertions). The first term on the RHS is the expectation value of the change in measure 
(with $k insertions) and the second term on the RHS is the expectation value of the change 
in insertions <&k ■ So far, we have not made any approximations. Let us now expand the multi- 
trace correlations according to and the factorization formula ([3]). The SDE at order 1/N° 
are the large- N factorized SDE(fSDE) or loop equations(LE). They only involve the large- N 
limits of single trace correlations Gj 

vfS M Gj lIj2 = v\b\^G h G h = vtfj. (8) 

Using the notation L v = v\L\ for the vector fields associated to the infinitesimal changes (Jl]), 
these LE may be written L v S J Gj = vfrfj where rfj = jj 1 * /2 G/ 1 G/ 2 . Here the action of the 



7 These changes of variable don't always preserve hermiticity of Ai . Under a change of integration variable, 
the integrand, measure as well as domain of integration may change, but the value of the integral is unchanged. 
It is possible to derive the SDE by making hermitian changes of variable, see appendix [XI 

8 {&k i • • ■ &K n ) is the quotient of two integrals. Here we make a change of variable in the numerator but not 
the integral for Z in the denominator. We could change variables in each, but this doesn't give new equations. 

9 For the purpose of deriving the loop equations and Ward identities, it is adequate to start with SDE for single 
trace correlations, so the reader could set n = 1 in a first reading. 
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vector fields on the moments is L\Gj = 5 J j llj2 Gj 1 u 2 and extends by linearity and the Leibnitz 
rule to polynomials in the Gj . Moreover, the Lie bracket of two such vector fields is 

[L},4] = S J /^L^ Ij2 -S I I ^L} lJl2 (9) 
or [L U ,L V ] = L w where L w = wj? L k K and 

<= E («i v k liK2 - u f liK2 <)■ (io) 

K=K 1 IK 2 

At 0(1/N 2 ), the SDE (one for each v and n > 0) involve the Gj as well as the G^ 's 

„,I aJiiJ2r<( 2 } — r,J xhihM 2 *) 

n 

+ v i E PG ^i ' ' ' g k p „ 1 Gl p im p Gk p+1 ■ ■ ■ G Kn . (11) 
P =i 

We could continue listing the SDE at each order in 1/-/V 2 , but we refrain from doing so since 
they no longer involve the single trace correlations Gj . Gj are the primary objects of interest in 
the large- N limit and our goal is to determine them for a given action S(A) . It would be ideal 
if we could uniquely determine them by solving the LE ([8]) . Unfortunately, as was demonstrated 
in [20] (and reviewed in !2.2p . this is not possible for many interesting actions S(A) , since the LE 
are underdetermined. One source of this problem was that there are vector fields v for which 
both LHS and RHS of ([8]) identically vanish for al0 Gj , so that the LE for those v are vacuous. 
Such vector fields are associated to simultaneous symmetries of the action and measure in the 
large- N limit. We will call such symmetries non- anomalous symmetries of the large- N limi10- 
Of course, in general, v need not be a symmetry of either action or measure. 

We would like to use the 0(1/N 2 ) SDE CD]) to determine the Gj that the LE do not fix. 
In principle, (jl If) are always valid. However, (llip involve the G^ 's which we do not wish to 
determine (and most likely cannot, without also involving the O(^) SDE and so on). Thus, we 
would like to use the subleading SDE (lllj) only for those v for which the LE (|8|) are vacuous. But 
even these equations would seem to involve the pesky G^ 's. Fortunately, a remarkable stroke 
of good fortune comes to our rescue. Suppose a vector field v is such that it is a simultaneous 
symmetry of the action and measure at N = oo, i.e. vj S JiiJ2 Gj 1 jj 2 = = v-rfj for all Gj. 
Then we will show (sections I3.3I and I4.1jl that v is also a simultaneous symmetry at finite iV , 
and thence a symmetry at each order in 1 /N 2 : 

vj6 I / iI *{<!> h <S> l2 $ Kl ---<f> Kn )= =v{S- hiJ H<S> JlI j 2 <5>K 1 ---<S>K n ) V n = 0,1,2,- •• 

if vlS J ^G JlIj2 = = vjr,} V Gj. (12) 

Thus, the terms involving the G^ 's in (|lip would identically vanish for such v and (lllj) would 
reduce to a set of Ward' identities 

n 

viY, 5L K p Mp GK 1 ---G Kp _ l G LpIMp G Kp+1 ---G Kn =0 V Kj, n > 1. (13) 
P =i 

10 When we say 'for all Gj ', we really mean 'for all cyclic and hermitian Gj '. 
11 Anomalous symmetries leave the action invariant but not the measure. 
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We call these 'Ward' identities (WI) since they are analogues of the Ward-Takahashi-Slavnov- 
Taylor identities of Yang-Mills theory. The latter are a consequence of BRST changes of variable 
(0) in functional integrals. Recall that the BRST transformations are also non-anomalous in the 
sense that they leave both the gauge fixed Yang- Mills action and measure invariant. 

These Ward identities can be written more compactly as 

L v {Gk 1 ■ ■ ■ Gk„) = V K p , n> 1. (14) 

Those for n > 1 follow from those for n = 1 and the Leibnitz rule. So the WI may be taken as 

L v G K = V K provided v is such that v\ S Jllj2 G JlIj2 = = v{rj\ V Gj. (15) 

It is satisfying to see that WI, which arise as a consequence of non-anomalous symmetries, may 
be regarded as a special case of the more general concept of Schwinger-Dyson equations. This 
is really a statement about quantum field theory in general, though we are discussing matrix 
models here. Traditionally [23] Ward-like identities are not regarded as related to Schwinger- 
Dyson equations in this manner. As pointed out in [20J and reviewed in section [221 the factorized 
large- N SDE are often insufficient to determine the correlations of a matrix model. However, 
when the fSDE are supplemented by the above WI, it becomes possible to determine many (and 
possibly all) the correlations, as we will see in later sections. 

2.1 The case of a single matrix 

For a 1 -matrix model with action tr S(A) = tr J2i<n<m S n A n , we use the changes of variable 
L v : A — > A + Xm>-i v n A n+l . A convenient basis is L n A = A n+1 ,n = —1, 0, 1, ... . These are 
familiar from the Lie algebra of polynomial vector fields on the real line, L n = x n+1 J|,re = 
— 1, 0, 1, Their Lie bracket is [L m , L n ] = (n 

L u = XI UnLn with i L u,L v ]=L w where w k = X (n-m)u m v n . (16) 

n> — 1 m-\-n=k 

m,n> — 1 

Their action on the moments G n = lim^^co <$„> is L k G n = nG k+n . Here $ n = %A n . The 
moments are real by hermiticity of A. The LE are 

m 

X] v k ^nS n G k+n = X v k 7 lk = X Vk X G p G q- ( 17 ) 

fe>_l J1=0 fc>-l fc> — 1 p+q = k 

p,q>o 

The LHS is the expectation value of change in action J2 k >-i v k,L k J2i<n<m S n G n while the RHS 
is the expectation value of the infinitesimal change in measure in the large- iV limit. 

The only vector fields for which the N = oo expectation value of the change in measure 
vanishes, are translations A — > A + v _il . To see this note that the change in measure term is 

Y, v k X G P G <i = v o + 2v x G x + v 2 (2G 2 + G\) + ■ ■ ■ (18) 

k> — 1 p+q=k 
P,9>0 

If this is to vanish for arbitrar\l^l real G\ , G% , ■ ■ ■ , then we must have vq = v\ = v 2 = ■ ■ ■ = ■ 
Thus only V—\ can be non-vanishing, which corresponds to a translation. 

12 The only constraints on the moments are that they be real and satisfy the moment inequalities, i.e. that the 
Hankel matrix gij — Gi+j be a positive matrix. 
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The only action for which translations are a symmetry in the large- N limit is the trivial 
action, S(A) = constant: the expectation value of change in action under a translation is 

m m 

v k Y nS n G n+k = v-\ Y nS n G k+n = v_i(Si + 2S 2 Gi + 3S 3 G 2 H h m5 m G m _i). (19) 

k>-l n=0 n=0 

If this must vanish V G n , we must have S\ = S2 = • ■ ■ = S m = 0, i.e. a trivial action. Thus, 
for a 1 -matrix model, we have no infinitesimal simultaneous symmetries of action and measure 
in the large- N limit. Consequently, there are no WI to supplement the LE with. 

This leaves a small mystery for 1 -matrix models. As discussed in [20], the LE (|17p of a 1- 
matrix model are underdetermined. They do not fix G%, G2, ■ ■ ■ G m -2 ■ The higher moments are 
fixed in terms of these by the LE. How are the first few moments to be determined if there are 
no WI to supplement the LE? Of course, for a 1 -matrix model, there are alternative techniques 
such as solving the integral equation for the eigenvalue density [2]. For multi-matrix models, 
the LE are often more severely underdetermined (there are an infinite number of moments that 
are not fixed). Remarkably, for multi-matrix models, where no alternative systematic method 
of solution exists, the WI do alleviate the underdeterminacy of the LE (section [5]) . 



2.2 Underdeterminacy of multi-matrix loop equations 

The multi-matrix LE (|8|), can also be written as S JllJ ' 2 Gj 1 u 2 = 5 I I lll2 Gi 1 Gj 2 for each i and 
/. This form, where we take the monomial basis L\ for vector fields L v is convenient for our 
current discussion. In general, these LE are underdetermined, as found in section 2.2 of Ref.|20j. 
Part of the reason for this underdeterminacy is the presence of non-anomalous symmetries of 
action and measure. First, we establish that the LE for given /, i can be regarded as a system 
of inhomogeneous linear equations for higher rank correlations with lower rank ones possibly 
appearing non-linear ly. From the LE, it is clear that if there are any correlations appearing 
on the LHS, they will be of a higher rank than the ones on the RHS. More precisely, suppose 
the action is an m th order polynomial (i.e. there is a non-vanishing coupling tensor S K with 
\K\ = m). LHS of the LE for given I,i (if it is non-trivial^!) , involves correlations only linearly 
and with a rank between \I\ and \I\ + m — 1, while the highest rank correlation on the RHS 
has rank \I\ — 1 . Even if the LHS vanishes, the highest rank correlation in the LE still appears 
linearly, but now on the RHS, and has rank \I\ — 1. However, in many cases, we find that this 
system of linear equations is inadequate to determine all Gi . 

Let us illustrate this with a Gaussian + Yang-Mills 2-matrix model tr S(A) = tr [y(vli + 
A2) — ^-[A\, A2} 2 ] ■ The matrix integrals for this model converge, and the correlations make 
rigorous sense and could for instance be measured numerically. So it makes sense to try to find 
them by solving the LE. In this case, the cyclically symmetric coupling tensors are 

5-H = S 22 = ^, S 1122 = S 2112 = S 2211 = S 1221 = l/(4a) and S 1212 = S 2121 = -l/(2a). (20) 



13 Even if there is an S K 7^ with \K\ = m, it may still happen that for some choice of I and i, the coefficients 
of all correlations of rank |/| + m — 1 on the LHS of the LE vanish. An example is the Gaussian+YM 2-matrix 
model LE with m = 4 and empty I , given later in this section. Thus, it is not true in general that the maximal 
rank correlation appearing in a LE has rank I\ + m — 1 . This possibility, which is special to multi- matrix models 
and has no analogue for 1-matrix models was overlooked in Ref. |2U] . 
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The LE are 



i = 1 : mG n - a-\2G I212 - G I221 - G I122 ) = 5j llh G h G l2 and 

i = 2: mG I2 -a- 1 (2G I1 2i-G I n2-G I211 ) = 5 I I l2l2 G h G l2 . (21) 



For 7 = 0, the LE sajo G\ = G 2 = 0. The LE for |/| = 1 relate 2- and 4-point correlations: 

1 = 1: mGn - -(2G m2 - 2G 1221 ) = 1 & mG 12 - -(2G U21 - G 1U2 - G 12U ) = 
a a 

1 = 2: mG 21 - -(2G 22l2 - G 2221 - G 2122 ) = & mG 22 - -(2G 1212 - 2G 1221 ) = 1.(22) 
a a 

They give the conditions Gi 2 = G 2 \ = 0, G\\ = G 22 and mG\\ = 1 + -^{2G\ 2 i 2 — 2G\\ 22 ) . They 
do not determine G\\ and give only one relation among the 6 independent rank-4 moments. 
LE with |/| = 2 relate 3- and 5-point correlations (since we already found Gi = .) 

7 = 11: mGin - -(2Gn2i2 -2G m22 ) = & mG n2 = 0. 
a 

I = 12 : mG\ 2 \ (G12122 — G\\ 222 ) = & mG\ 22 \G\\ 2 \ 2 — G11122) = 

a a 

I = 21 : mG 2 u (G12122 — ^11222) =0 & mG 2 \ 2 (G11212 — G11122) = 

a a 

1 = 22: mG 221 =0 & mG 222 - -{2G l2122 - 2G ll222 ) = 0. (23) 

a 

The |/| = 2 LE imply that all G^ vanish and give two relations among the 8 independent!^! 
5th rank moments, G\ 2 \ 22 = G\\ 222 and G11212 = Gni22- The |/| = 3 LE relate 2- 4- and 
6-point moments (we omit those equations that contain no new information) 

1 = 111: mG 1 iii--(2Gi 112 i2-2Gi 11 i22) = 2Gn and mG 1112 = 

a 

I = 112 : 2G112212 = G\i\ 222 + G112122 

"lGn22 (2Gni212 — G112112 — G111122) = Gu 

a 

I = 121 : 2G121212 = G112122 + G112212 

2 

"1^1212 = — (G112112 — G111212) 

a 

I = 211 : 2G112122 = G112212 + G111222 

1 = 222: mGi 222 = and mG 2222 (2G121222 - 2G112222) = 2G 2 2 

I a 

I = 122 : mGn22 (2G121222 — C112222 — G122122) = C22 

a 
2 

7 = 212: mGi2i2 = — (C122122 — ^121222)- (24) 

a 

However, these 11 equations (even if all are independent), are inadequate to find the c(n = 6, A = 
2) = 14 independent rank- 6 correlations, let alone the unknown 2 & 4-point correlations. 

Similarly, consider a Yang-Mills 2-matrix model tr S(A) = — ^— tr L4i,^2] 2 - The matrix 
integrals do not converge here due to the flat directions in the commutator squared action. So 



14 This example illustrates that even in a quartic model (m=4), the LE may determine correlations of rank 
m — 2 = 2 or less. 

c(n, A) denotes the dimension of the space of cyclically symmetric hermitian tensors of rank n in a A matrix 
model. A formula for c(n, A) is given in appendix iDl 
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our derivation of the LE and WI are not strictly valid in the case, though they can be considered 
formally. In particular, it is not clear that the LE and WI form a consistent system of equations 
for this model. Nor is it clear how one could check an answer for a particular correlation, say 
by Monte Carlo integration, since the matrix integrals do not converge. Nevertheless, we can 
consider the LE formally here in order to show that they are under determined. The LE are 

-vf (G/122 + G/221 - G/212) + -viiGnn + G m2 - G I121 ) = [v[5^ lh + vffi^lGj.G^. (25) 
a a 

Since vj are arbitrary, we get a pair of LE (for each word / with |/| > 0) 

G/221 + G/122 - 2G/212 = ad 1 / 112 G ^Giz, G/112 + G7211 - 2G/121 = ad^^G^G^. (26) 

All correlations of rank 1 or 2 are undetermined. In addition, taking 7 = does not give any 
relation for third rank moments, since the LHS of the LE identically vanish on account of cyclic 
symmetry. As for rank- 4 moments, we get only one relation 2G1212 — 2G1122 = —ct, from the 
LE, which is inadequate to fix the 6 independent 4 th rank correlations. 

Similarly, the LE of the Chern-Simons 3-matrix model tr S(A) = 2«g*i fc tr AiAjA^ and 
Mehta 2-matrix model [M] tr \cA x A 2 + (g/A)(Af + A$)] are underdetermined (§2.2 of [20]). 



3 Measure preserving transformations 

Our aim in this section is to determine the vector fields L v : Ai — > Ai + v\ Aj under whose 
action the matrix model measure is invariant. We call such transformations measure or volume 
preserving. These vector fields are universal in the sense that they are independent of the choice 
of action S(A) . They can only depend on the size of the matrices ( N), the number of matrices 
(A) and on the ensemble from which the matrices are drawn (hermitian in our case). 

The main result of this sectional is that vector fields L v = v\ L\ satisfying (|55p are measure 
preserving for any N . In the large- N limit these are the only ones, but for finite N there could 
be more. In particular, a symmetry of the matrix model measure at N = 00 is automatically 
a symmetry of the measure at finite N and consequently at each order in 1/N 2 . A simpler 
sufficient condition for a vector field to be measure preserving is 

v p ( j) + v (j W) = 0) v x j (27) 

Here (• • •) denotes cyclic symmetrization ()38|) . Measure preserving vector fields form an infinite 
dimensional Lie algebra for A > 1 with Lie bracket © (see section I3.4p . For A = 1 , it is a 
1 -dimensional abelian Lie algebra consisting of translations A — > A + (see section l2"TTj) . 

3.1 Change in measure due to action of (homogeneous) vector fields 

If all I appearing in the components v\ of the vector field L v = v\L\ have the same length 
\I\ , then we will call such a vector field homogenous of rank |/| . The variation of the measure 
under (a not necessarily homogenous) infinitesimal change of variable, 

A i ^A' i = A i + v\A I (28) 

16 Section l3.2l where this is established is a bit long and can be skipped in a first reading. 
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is the first order term in the expansion of the determinant of the Jacobian J in powers of v\ 
dA't 



J = det 



Here $>/ 



det^-Ot + w i J <y? ifa *i 1 $j a + 0{v z )} = 1 + N z vl5Y ll2 <5> h <$> l2 + 0{v z ). (29) 



2„J X hih; 



SL A T 



Thus infinitesimally, the change in the measure per N 2 is 



SJ_ 

N 2 



(30) 



We want to determine those v\ 's for which 5J/N 2 = for all <&k which are cyclically symmetric 
in K and hermitian <§>* K = $>^. So we should set the coefficients of independent $k to 
zero. Unfortunately, for finite- JV, the analysis is complicated by the fact that are not all 
independent. Indeed, they are related by trace identities (analogues of conditions from vanishing 
of characteristic polynomial for a single matrix). However, in the large- N limit, the #/ (and 
consequently their expectation values, Gj = (<J>j) ) are independent up to cyclic symmetry and 
hermiticity. The expectation value of the infinitesimal change in measure becomes 



lim (-^) 

N^oo^N 2 



Gii Gi 2 



v*4 



v ■ rj. 



(31) 



Setting the coefficient of each independent monomial in the Gi to zero leads to a characterization 
of the measure preserving vector fields v\ in the large- N limit (section I3.2[) . 

An inspection of v ■ rj reveals that it vanishes under the sum of two homogeneous transfor- 
mations 5A{ = J2\i\=const v i Ai + J2\ j\=const v i Aj with |/| ^ \ J\ if and only if it vanishes under 
each separately. So without loss of generality, we restrict to homogenous vector fields. 

Notice from ([9]) that the commutator of homogenous vector fields of rank p > 1 and q > 1 is 
a homogenous vector field of rank p + q — 1 . Using this, we define the grading of a homogenous 
vector field Y^\i\=const v i^i as l-^l — !■ With this, the Lie algebra C of all L ' J^l becomes a Lie 
algebra C = p >o £p graded by the non- negative integers, where 



C p = span |l} : |/| = p + 1 j and [C p , C g ] C C p+g . 



(32) 



This can be used to generate homogeneous higher rank volume preserving vector fields from 
ones with lower rank, provided the latter do not form an abelian Lie algebra. 



3.2 Characterization of measure preserving vector fields for N = oo 

Roughly, the condition that a vector field be volume preserving becomes stronger as N — > oo, 
since the number of independent trace invariants grows in this limit. So there are potentially a 
lot more volume preserving vector fields at finite- ./V than at N = oo. Fortunately, the condition 
that a vector field be volume preserving at N = oo will be seen to be a sufficient condition for 
it to be volume preserving at any finite N . In this manner, we will establish that if the change 
of measure term in the LE ([8]) vanishes for a given vector field v , then it also vanishes^! in the 
finite- N Schwinger-Dyson equations ([7]) and indeed at each order in 1/N 2 for the same vector 

17 Except translations v\L\ : Ai — > Ai + vf 1 , which are a separate abelian algebra. is the empty word. 
18 There likely exist vector fields preserving the measure at finite- TV but not at N — oo. 
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field v . To characterize volume preserving vector fields in the large- ./V limit we must solve the 
equations v ■ rj = for v\ . Let us begin with homogeneous vector fields of lowest rank. 

• Constant shift: SAi =v{S.. In this case v ■ rj = 0. So all homogeneous vector fields of rank 
zero are symmetries of the measure. This reflects translation invariance of the measure. 

• Linear transformation: SAi = v\Aj are measure preserving iff they are traceless: 

v ■ rj = vfSj = v\ so !)'j} = <S> d is traceless v\ = 0. (33) 

• Quadratic: SAi = AjA^ . In this case, 

v V = v i\6lG J +8p k )=vrG j +vl k G k = ^ + vr)G J (34) 

Thus quadratic vector fields that preserve the measure must satisfy v % ? + v\ l = . 

• Cubic: SAi = v{ kl AjA k Ai 

This must vanish for all cyclic and hermitian Gk ■ The linear term in G 's is cyclically sym- 
metric in mn, so it is not necessary that the coefficient of G mn and G nm separately vanish. 
Rather, only the cyclic projection of G mn 's coefficient must vanish. Similarly, the quadratic 
term in G's is symmetric under m^nso only the symmetric projection of its coefficient must 
vanish. Moreover, hermiticity implies G* = Gj and G* mn = G nm = G mn , so all 1- and 2-point 
correlations are real. We need not worry about setting the coefficients of their imaginary parts 
to zero. Thus v ■ r\ vanishes identically if and only if 

v nnn +v mrn +v mni +v nmi = Q and v min + y nim = Q ^ 
We can write this more succinctly as 

v Kmn) +v (mn)i = Q ^ v min + v nim = Q (37) 

Here we introduced the cyclic symmetrization operation (• • •) which is defined aJ^l 

U)iK \ - n(J)iK Ji(K) \ - Jia(K) (J)i(K) \ - nU)ia(K) , 00 x 

7reqj| creq K | ^ c \j\ 

<reC|K| 

jklm 



Quartic: SAi = vj AjA k AiA r/ 



V 'V - V l m ^Jklm G Pqr + SP jkT m G P G qr + ^jklm G P1 G r + ^jklrrPvqr) 

= (v7 qr + v™ H )G pqr + vf qr G p G qr + vr r G pq G r 

= w qr +v!n G p q r+(vf iqr +vrnG p G qr . m 

Here, the quadratic term is cyclically symmetric in qr, so its coefficient must be cyclically 
symmetrized in qr . Similarly, the linear term is cyclically symmetric in pqr so we must cyclically 
symmetrize its coefficient in pqr. However, there is a further subtlety that we must address: 



19 C|j| is the cyclic group of order J . Note that we do not divide by the number of terms. 
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Gk are complex numbers, but their real and imaginary parts are related to those of via 
G* K = G^ . Hermiticity implies that G p and G qr are real, so for the quadratic term, it is 
necessary and sufficient that (vf are real) 



v pi{qr) _|_ y (qr)ip 



0. 



On the other hand, G* 



Grqp • So ^R.Gpg r 



ffiG^qp and ^sG pq r 



(40) 

QG rqp . So it is necessary 



and sufficient to set the coefficients of $lG pqr + $tG rqp and ^sG pqr — QG rqp to zero separately: 



»(P9T) i (P1 r )i i i( r QP) i (»"9P)» 



n J *(P<3 r ) i ifrgp) (''9P)j 

and Va + t>- — v A — v) 



V pqr. (41) 

However, these two conditions are equivalent (by adding and subtracting) to the single condition 

v fPV) +v tor)i = Q V pqr. (42) 
Therefore v ■ r/ vanishes if and only if 

4 pqr) + vl pqr) * = and vf qr) +v ( - qr)ip = 0. (43) 

• Quintic: For a rank 5 vector field A{ — > j4j + uf <?rs *^4p gr st to be volume preserving we need 



V ■ T) 



pqrs 



ipqrs 



+ v pqrsi )G pqrs + {v piqrs + v qrsip )G P G qrs + {v pqtrs + vl sipq )G pq G rs = 0. (44) 



Since moments of different ranks are independent, v ■ rj = iff the following three equations are 
satisfied (we have cyclically symmetrized as in previous cases in order to reduce to a sum over 
equivalence classes under cyclic symmetry, which is denoted ~) 



E <<■■■ 

pqrs/r- 



')('j„ir, y2 

p,qrs/~ 

and 2_, 

pq/~,rs/ 



v pi{qrs) _|_ v {qrs)ip 



, (pq)i(rs) 



)G p G qrs 

(rs)i(pq) 



0. 

)G pq G rs 



(45) 



It remains to take care of the relations imposed by hermiticity to select the independent mono- 
mials. Consider the first equation in (|45p. By hermiticity ?R.G pqrs = $lG srpq and QG pqrs = 
—QG srqp . So we further restrict the sum to equivalence classes under reversal of order of in- 
dices. We will denote the combination of the quotient by cyclic symmetrization and reversal of 
order of indices by the symbol ~' . Then the first condition in (I45j) becomes the pair 



E 

pqrs/^j' 

E 

pqrs/^' 



i(pqrs) (pqrs)i i(pqrs) 



(pqrs)i 



RG 



pqrs 



y i(pqrs) _|_ v (pqrs)i _ ^(pqrs) _ v (pqrs)i 



%G 



pqrs 



and 
0. 



Since the sum is over independent moments, we set the coefficients to zero and get 



v 



i(pqrs) 
i 

i(pqrs) 



(pqrs)i . i(pqrs) 



(pqrs)i 



% ' i 1 i 



and 
0, 



(46) 



(47) 
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for each equivalence class pqrs under the relation ~'. However, this pair is equivalent to 

V cyclic equivalence classes pqrs/ ~ . (48) 



v i(pqrs) _|_ v (pqrs)i 



A similar analysis of the last two conditions in (|45p using hermiticity (G* = G p and $tG qrs = 
WtGqrs and QG qrs = —^sGqrj and G* q = G pq ) allows us to identify the coefficients of independent 
moments and set them to zero. When the dust settles, the necessary and sufficient conditions 
for a 5 th rank tensor to be measure preserving are 



i(pqrs) , (pqrs)i 
v i + v i 



0; 



v pi(qrs) _|_ v (qrs)ip 



0; 



v (pq)i{rs) _|_ v (rs)i(pq) 



Vp,q,r,s. 



(49) 



We see that so far, the hermiticity relations between the Gk, though taken into account, did 
not make their presence felt in the final answer. This simplification is due to Gi and Gij being 
real. The hermiticity relations will play a role in the necessary and sufficient conditions for rank 
7 and higher vector fields to be measure preserving. This is because it is the first case where r\\ 
involves quadratic monomials in moments where both factors can be complex, e.g. G pqr G s t u - 
This leads to complications which we now deal with in the general case. 



Rank n : In the general case, 8A{ 



Aj 1 . . . Aj n and v ■ rj 



Mill 



Gi 1 Gi 2 is a quadratic 



polynomial in moments. The necessary and sufficient conditions on v for v ■ r] = are got 
by selecting the independent monomials and setting their coefficients to zero. We use three 
relations: (a) commutativity of products GiGj = GjGi , (b) cyclicity Gi = Gj if I, J are 
cyclically related and (c) hermiticity Gi = G*j or $lGi = KG/ and QGi = —QGj. (a) we 
must symmetrize the coefficients in I\ and I2 and restrict the sum over I\ and I2 to include 
only (say) Gi 1 Gi 2 and not Gi 2 Gi 1 (this is denoted X) )■ 



v ■ rj 



& 



Y.\v{ lll2 +v[» h )G h G l2 
h,h 



0. 



(50) 



Relation (6) means we must cyclically symmetrize coefficients in I\ and I2 and further restrict 
the sum to cyclic equivalence classes of I\ and I2 (denoted 1%/ ~ ) 



v ■ r\ 



* 



h/- 



v (h)i(h) +v (h)i(h) 



Gj 1 Gj 2 = 0. 



(51) 



Implementing (c) is more tricky. We must identify monomials that are independent after ac- 
counting for hermiticity. Taking !R & 9 parts (v( € R), we write v ■ r) = as the pair 



E 
E' 



(KG/KGj - 9G/9Gj) 
(sRG/SjGj + KGjQG/) 



and 



(52) 



The last two terms can be combined. Hermiticity = QGi = —QGj. So tftGiRGj 

is independent of ^GjQGj and we can set each part to zero separately. Thus v ■ rj = iff 



E' 



..WW 4. w (-o*w 

i ' i 



31GjKGj = 



E' 



SGftGj = 0, 
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and ^ 



I/r. 

J/r 



m^Gj = 0. (53) 



Here the sums include words / as well as their mirror images /, so the monomials such as 
?R.Gj^R.Gj are not all independent on account of the hermiticity relations. We must further 
restrict the sums to equivalence classes under reversal of order of letters in a word to get a truly 
independent basis for quadratic polynomials. Once this is done we set the coefficients to zero 
and find that v ■ r] = iff (the signs are determined by the hermiticity relations) 

(i)i(J) . (J)i(J) . (i)i{J) . (i)i(J) . j j n 

(I)i(J) (I)i(J) (I)i(J) . (I)i(J) . j j n A 

v\ — v\ — v\ + v\ + 1 ^ J = and 

(i)i(J) . (J)i(J) (i)i(J) (I)i(J) . T T n fKA ^ 
v i + v l — v\ ' y ' — v\ + 1 <-> J = 0. (54) 

These can be slightly simplified to the following three conditions 

v\ +v\ ' + 1 <-> J = 0, v\ +v\ + I <-> J = 0, v\ v l + 1 J = 0- (55) 

Thus, a homogeneous vector field v\ of rank n is volume preserving at N = oo (v ■ r] = 0), iff 
conditions (|55p are satisfied for each multi- index / and J such that |/| + |J| = n — 1. Since 
conditions (|55p are somewhat lengthy (though easy to remember), it is pertinent to add that a 
sufficient (but in general not necessary) condition for v to be measure preserving is 

v (.mJ) + V (JW) _ o v /, J with \I\ + \J\=n-l. (56) 

More explicitly, this sufficient condition may be written as ( [n] is the greatest integer part of n) 

*0'l— Jn-l) I Ul-jn-l)i _ n 

u i ' u i u ' 

jli(h---jn-l) . (j2-j„-l)ijl _ n 

u i ' u i U J 

,,(jih)i(ja~'3n-i) , 0'3-J'n-i)»(3iia) _ n 

O'lia— irn-l])»0',n+lj—Jn-l) (j'r n+1 , — 3n-l)i(jl32—j,n-l ,) 

^ +Vi -a- = 0. (57) 

In fact, (|56|) is both necessary and sufficient for vector fields of rank < 6. 

Now, it is easy to see that a volume preserving vector field for N = oo is automatically volume 
preserving for finite N. Suppose v\ is such that v-r\ = v\Sj l%l2 Gi x Gi^ = 0, or equivalently, such 
that conditions ([55]) are satisfied. Then it will automatically satisfy Mx = vjSj 1 ll2 $7-^/2 = 0. 
For, all we needed for v ■ r] = was commutativity of the product of two G\ 's, cyclic symmetry 
and hermiticity of the Gi . All these properties are also satisfied by the 's. Of course, there 
are likely to be vector fields other than those satisfying (I55p (i.e. v ■ 77 ^ 0) for which jfc = 0. 



3.3 Volume preserving vector fields annihilate measure terms in LE and SDE 

Thus far, we have shown that vector fields characterized in (|55h leave the measure invariant 
5 J/N 2 = v I i 5 I j 1%l2 ^i 1 ^i 2 = both for finite and infinite N. Moreover, they were all the vector 
fields that left the measure invariant for N = 00: v{5 I I lll2 Gj 1 Gi 2 = vjrji = 0. In other words, 
the RHS of the LE ([8]) identically vanish iff the vector field v satisfies (|55p . 
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On the other hand, multiplying ([30]) by $>Ki • • • $K n an d taking expectation values we get 

vffi Ua {$ Il $ Ia $ Kl ---$ Kn ) = Q V $ A - and n = 0,1,2,... (58) 

provided v satisfy (|55|) . Combining with the result of the previous paragraph, we see that vector 
fields for which the RHS of the LE ([8]) vanish, also annihilate the change of measure term on 
the RHS of the finite- iV Schwinger-Dyson equations ©. Furthermore, multiplying by iV 2 and 
letting N — > oo we see that the same class of vector fields also annihilate the change of measure 
term on the RHS of the 0{l/N 2 ) SDE ^ 

vW U2 G% h]Kl ,., Kn = 0. (59) 

This is a part of the result we needed in section [2] to establish the WI (|15p . The other part 
involves identifying which of these volume preserving vector fields also leaves the action of a 
specific matrix model invariant, a task we will undertake in section 

3.4 Volume preserving vector fields form an infinite dimensional Lie algebra 

It should be possible, but laborious, to check that the Lie bracket of two vector fields of the 
form (|55p is again of the same form (we have checked this for vector fields of some low ranks). 
But there is a simpler argument (which uses a much deeper result from |13j) that shows they 
form a Lie algebra. In [13J it was shown that there is an entropy functional x such that 
L\x = rfi = 8j 1%Ia GiiGi2- Suppose L U ,L V are volume preserving. From our results in sections 
I3.2l and l3.3l this means u\r\\ = and v^rjj = 0. Then we have L u \ = and L v \ = 0. It follows 
therefore that [L u ,L v ]x = (L U L V — L v L u )x = 0. Thus L w = [L U ,L V ] is also volume preserving. 
We conclude that volume preserving vector fields form a Lie algebra. 

Example: The measure preserving vector fields corresponding to linear transformations, 
L u = UjL? form the sZa(R-) Lie algebra for a A-matrix model. We already found (section 13. 2p 
that measure preserving linear transformations are the traceless ones. Here, we check that their 
Lie bracket implied by Q is the same as the sZa(R-) Lie algebra. 

[L U ,L V ] = u}vf[Li,L{] = u)vl(5{L\-5\L{) = u'^-ujw*^ = w\L\, (60) 

where w\ = uj.vf — v\v§ = {[u,v\)\ is just the matrix commutator. Thus [L U ,L V ] = L\ u ^] and 
the linear symmetries form the Lie algebra sl\(R) . 

Moreover, for A > 1 we can show that the space of measure preserving vector fields is infinite 
dimensional. It is sufficient to consider each rank separately. First, the space of rank-n vector 
fields y l . lt2 '" ln is A n+1 dimensional. For a rank-n vector field to be measure preserving it is 
sufficient (though not necessary) that it satisfy equations (f56j) . There are at most A™" 1 such 
linear equations (if they were not linearly independent or necessary, there would be even fewer). 
Thus, the space of solutions is at least A ra+1 — A n_1 dimensional. This grows exponentially with 
rank, so measure preserving vector fields are an infinite dimensional Lie algebra for A > 1 . 

20 However, \ cannot be expressed as a formal power series in Gi . 
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3.5 Explicit examples for 2 and 3 matrix models 



From sections 13.21 and 13.41 we know that linear volume preserving vector fields are traceless 
matrices vj , i.e. elements of sl\(R). This is a A 2 — 1 dimensional space (3 dimensional for a 
2 -matrix model and 8 dimensional for a 3 -matrix model). 

A generic quadratic vector field v{ k in a A-matrix model is specified by A 3 parameters. But 
volume preserving vector fields obey relations given in section 13.2^ which restrict the number of 
independent coefficients. Let us work out volume preserving v{ k for 2 and 3-matrix models and 
determine the dimension of the space of such vector fields. The condition for v\ to be volume 
preserving is + vj l = 0. In a 2 -matrix model this is the pair of equations 



2v\ x + vf + vf 



and 



2vf 



.12 



+ V 



21 



0. 



(61) 



So quadratic volume preserving vector fields are the 2 3 - 2 



L, 



„,llr2 , 7 ,22rl ,12 
v 2 -^11 T V 1 L<22 T V 1 



+ 



vf 



L 



12 



-L 



ii 



+ vf 



12 

-2 
J 21 



-L 



6 parameter family 
1 



22 



^21 



-L 



22 



-L 



ii 



(62) 



In a 3-matrix model there are three independent conditions 



2v\ l + vf + + v\ 2 + vl 3 



0.. 



2vf + + vf + < + vi 



.12 



32 



.21 



.23 



and 2vf + + vf + vf + vf = 0. 
So quadratic volume preserving vector fields are a 3 3 



(63) 



3 = 24 parameter family for A = 3. 



4 Transformations that also preserve action 
4.1 Establishing validity of Ward identities: last step 

So far, we have identified the vector fields L v = v\L\ which leave the measure invariant in the 
large- N limit and observed that they continue to be measure preserving even at finite N . In 
order to obtain the WI (|15j) , we need to determine which among these L v are also symmetries of 
the action. These are the non-anomalous infinitesimal symmetries. The answer will, of course, 
depend on the action of the matrix model being studied. For the infinitesimal change in ( x ) 
the action S(A) = tr S J Aj to vanish under A, — > Ai + v(Aj, we need 

L v S J §j = vl S' hiJ ^ JlIj2 = 0. (64) 

However, for finite N , not all the are independent even after accounting for cyclicity and 
hermiticity, due to the trace identities and other such constraints satisfied by the . So it is 
not straightforward to identify the necessary conditions on v\ . But in the large- N limit we 
may treat the as independent variables (up to cyclicity and hermiticity) . Taking expectation 
values, we must solve for v\ in the equations 

v!S Jlij2 G JlIj2 =Q V cyclic and hermitian G/. (65) 
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For such vector fields, the LHS of the LE ([8]) identically vanish. Moreover, a vector field that 
solves (|65p ) will automatically solve the finite- N equation (|64p ). though the converse need not 
be true. This is because all we use is cyclicity and hermiticity of Gj , which is also true of the 
<]?/ . Now multiplying (f64"P by $x x • • • &k„ , the same vector fields also satisfy 

vjS^Qj.u^K,---®^ =0. (66) 

Taking expectation values, we see that symmetries of the action in the large- N limit automati- 
cally annihilate the change in action term (with insertions) appearing in the finite- N SDE ([7]). 
In particular, multiplying by iV 2 and letting N — > oo, we see that the vector fields satisfying 
(|65|) also annihilate the change of action term on the LHS of the 0(1/N 2 ) SDE (|lip 

v [S J ^G% 2Ki ... Kn = 0. (67) 

Combining this with our result from section (|3.3p on volume preserving vector fields, we 
come to the following conclusion. Suppose the vector field v is such that both LHS and RHS of 
the large- N LE identically vanish, v I i S Jllj2 Gj 1 jj 2 = = vfrjj. Then the change in action and 
change in measure term in the 0{1/N 2 ) SDE also vanish identically v{ S^^Gj^j^...^ = = 
vfSf-^Gj ■i r K 1 —--K n ■ As a consequence, for such vector fields (non-anomalous vector fields), 
the 0(1/ N 2 ) SDE become WI (|15p which may be summarized as L v Gk = for all K . This 
completes the proof of validity of the WI. 

4.2 Non-anomalous symmetries of specific models 

It is straightforward to see that non-anomalous vector fields form a Lie sub-algebra of the infinite 
dimensional Lie algebra of measure preserving vector fields (section l3.4p . For, L V S(G) = and 
L W S(G) = implies that [L v , L W ]S(G) = 0. However, this Lie algebra is not necessarily infinite 
dimensional and depends on the action of the matrix model. 

This brings us to the task of determining the non-anomalous infinitesimal symmetries of 
specific matrix models. In looking for measure preserving vector fields, recall (section 13. ip that 
we could break up the problem into finding homogeneous measure preserving vector fields of 
a given ranlj 2 ^!. The same strategy does not work in general for symmetries of the action. 
However, if the action is itself a homogeneous polynomial 2 ^, then (|65p does not mix vector fields 
of different ranks. In that case, every solution to (|65p is a sum of homogeneous solutions. More 
generally, the action may not be a homogeneous polynomial, as for a Gaussian + Yang-Mills 
model. In such cases, not every solution of (|65p is necessarily a sum of homogeneous solutions, 
though there may still be large classes of homogeneous solutions. For this reason, we begin by 
determining non-anomalous homogeneous action-preserving vector fields of low rank. 

A priori, it is not clear that there are any vector fields that leave both action and mea- 
sure invariant. Indeed, for a 1-matrix model (section 12. lh there are none. We were pleasantly 
surprised to find not just linear but also non-linear non-anomalous symmetries for several inter- 
esting multi-matrix models. We begin with linear non-anomalous symmetries in section l4~3l and 
give examples of non-linear non-anomalous symmetries in section 14.41 

21 A homogeneous vector field v of rank n is one whose components vj are non- vanishing only for |J| = n. We 
call rank-1 vector fields linear transformations, rank- 2 vector fields quadratic changes of variable and so on. 
22 Examples include the Gaussian, Chern-Simons and Yang-Mills models. 
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4.3 Examples of linear non-anomalous symmetries 

We determine linear symmetries of both action and measure for the Gaussian A -matrix model, 
Chern-Simons 3-matrix model, Yang-Mills and Gaussian+YM A-matrix models. The linear 
non-anomalous symmetries of the Gaussian, CS 3-matrix model, YM 2-matrix model and Gaus- 
sian+YM 2-matrix models form the orthogonal Lie algebra with respect to the covariance matrix, 
s^R), sZ2(R) and o(2) Lie algebras respectively. Not every multi-matrix model has non-trivial 
linear non-anomalous symmetries. The pure-quartic 2-matrix model tr S(A) = tr (A 4 + i? 4 ) 
or the model studied by MehtajMj, tr S(A) = tr [cAiA 2 + (g/A)(Af + Af)] have no non- 
trivial linear action preserving symmetries. Linear symmetries form a closed Lie algebra among 
themselves (section [3.4p . so their Lie brackets cannot be used to generate new symmetries. 

4.3.1 Linear symmetries of Gaussian 

The Gaussian A-matrix model is defined by the action tr S(A) = tr ^C lJ AiAj where C %d 
is a positive real symmetric 'covariance' matrix. We seek infinitesimal linear transformations 
5Ai = v\Aj that leave the action as well as measure invariant in the large N limit. In section 
13.21 we found that the measure preserving transformations are the traceless ones v\ = forming 
the Lie algebra sIa(R) . Here we find that the vector fields v\ that preserve both the action and 
measure in the large- N limit are those that satisfy v l m + v k C kl Cj m = 0. This is the condition 
that v\ be an orthogonal transformation with respect to a metric given by the covariance. In 
particular, for a unit covariance C %i = 5 11 , these are the antisymmetric matrices. 

For the expectation value of the Gaussian action to be invariant at N = 00 , we need 

L V S(G) = v{C kl Gij = + v lj )Gij = V cyclic and hermitian Gy. (68) 

We used C kl and its inverse C; m to raise and lower indices, iP = v 3 k C kl , v jl Ci m = v 3 m , C kl Ci m = 
5™j . The condition for a symmetry of the action is that v be anti-symmetric 

xP l + v lj = 0. (69) 

If v is anti-symmetric, then v is automatically traceless. So action preserving linear transfor- 
mations are automatically measure preserving. To see this we first rewrite antisymmetry of 
v as a condition on v by lowering an index v l m + v k 'C kl Cj rn = 0. Taking the trace we get 
v\ + v 3 k C kl Cji = which implies v\ = 0. Thus, the non-anomalous linear symmetries of the 
Gaussian are given by vector fields v\ L*- that are anti-symmetric after raising an index with C %k . 
In other words, the orthogonal Lie algebra with respect to the metric given by the covariance 
matrix. In particular, the dimension of the space of linear symmetries ^(A 2 — A), does not 
change as we move around in the space of non-singular symmetric covariance matrices. 

Example 1: Consider a Gaussian two matrix model with diagonal covariance C 13 = diag(a, b) . 
Then the condition that 5Ai = vjAj be action preserving in the large- N limit is v \ = v\ = 
and av\ + bv\ = . These are the infinitesimal orthogonal transformations with respect to the 
'metric' diag(a, b) . Such v\ are traceless and thus measure preserving as well. 

Example 2: If the covariance of a Gaussian A-matrix model is a multiple of the identity, then 
the action preserving transformations are the anti-symmetric matrices (vj + v l j = 0) which form 
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the orthogonal Lie algebra with respect to the metric 5 13 . Such matrices are clearly traceless so 
that the non-anomalous linear symmetries form the Lie algebra o(A) . 

4.3.2 Linear symmetries of Chern Simons model 

The Chern-Simons 3-matrix model has action 

tr S(A) = —e^AAjAk = 2in tr Ai[A 2 , 4»] = 2ik tr (A 123 - A 132 ). (70) 

So the coupling tensors are S l ^ k = ^e l i k . "We seek all vector fields L v = v\L\ for which the 
change in ( jt x ) the action vanishes for N = oo (second equality requires relabeling of indices) 

L v S J Gj = v\S 3 ^G Jllj2 = 2iKv\^ k G l3k = 0. (71) 

vj are real and have no symmetry in /. They must satisfy vle l i k Gijk = for all cyclic and 
hermitian Gk- Specializing to linear transformations Ai — > Ai + vjAj, they must satisfy 

£ v\e^ k G ljk = 0. (72) 

l<i,j,k,l<3 

We could also arrive at this condition by making a linear change of variables in the action 

5-^trS = ^(viA i23 + v\A m + v\A m - v\A i32 - v\A m - v\A m ) = 0. (73) 
Writing out all the terms and using cyclicity of Gk this condition simplifies dramatically to 

3 

$>j (Cm - G 132 ) = 0. (74) 

8=1 

Taking real and imaginary parts^l we get the single condition v\ 9 G\ 23 = 0, which must be 
satisfied for all 3 G123 • We conclude that vj preserves the CS action iff it is traceless v\ = . We 
recall (section l3.2j) that traceless linear transformations also preserve the matrix model measure. 
Thus, the CS model has a maximal family of linear non-anomalous symmetries. 

From section E31 we know that the space of traceless real v\ is the Lie algebra ^(R), an 
8 dimensional space. The free parameters can be chosen as v\,v\,v\,v\,v\,v\,v\ and v 2 with 
v 3 = — v\ — v 2 . The corresponding symmetries are an 8-parameter family of vector fields 

L v = v\L\ + v\L\ + v\L\ + v\L\ + v\l\ + v\h\ + v\l\ + v\L\ - {v\ + v\)L%. (75) 

4.3.3 Linear Symmetries of Yang-Mills model 

For 2 or more matrices and a real symmetric invertible metric gij , the YM model has action 
tr S(A) = — tr [Ai, Aj][A k , Ai]g lk gi l . The expectation value of the change in the action under 
a linear transformation 5 A j = v\Aj in the large- N limit can be written as 

T Q(r<\ r< (n, m J-k„jl „,k iljm J im Jd , „,k Am (<7a\ 

L v b(G) - -—Gjkiml v { g g J -v { g g J -v i g g ■\-v i g g J J (7b) 
23 v\ € R. Hermiticity & cyclicity =>■ G123 = G132 which =>• KG123 = KG132 and QG123 = —QG132 . 
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To identify symmetries of the action, we must select independent Gjkim and set their coefficients 
to zero. First we restrict the sum to words jklm up to cyclic symmetry. Thus (55* = iff 

- £ G jklm Ri klm = (77) 



^ jklm I eye 



where the cyclically symmetric tensor R? klm is 



kirn 



vl(2g d g km - g im g kl - g lk g lm ) + cyclic(j - k - I - m -> j] 



{v {jl) g km + v {km) g jl ) - - 



yii m )g kl -|- y(3 k )g lm _|_ y( kl )g m J _|_ V ^ m ) gO k 



(78) 



Here we have used the metric to raise and lower indices v 3 = v[g and v 3 m = v 3l gi m and denoted 
the symmetric projection by v^ k ^ = ^(v 3k + v k3 ) . We have still to account for the hermiticity 
relations $tG jk im = ^Gjj^, QG jk i m = -^Gj^. Now, v\ and g kl are real and ^.G jk i m and 
QGj k i m are independent of each other, so we have L V S{G) = iff 

£ RP klm m jklm = and £ RP klm ZG jklm = 0. (79) 

jklm/cyc jklra/cyc 

Now we must collect the coefficients of 3tGj k i m and ^-G - klm and similarly for the imaginary 
parts and restrict the sum to avoid jklm if jklm has already appeared. Two possibilities arise: 
either jklm may be obtained from jklm via cyclic permutations or not. In the former case, 
jklm vanishes and the coefficient of §tGj k i m must vanish for v to be a symmetry of the action. 
Thus we get R? klm = if jklm is cyclically related to jklm. On the other hand, if jklm is not 
cyclically related to jklm, then collecting coefficients we have 



£ m 3klm {R ]klm + R jklm ) = and £ QG jklm (R 3klm - W klm ) = 0(80) 

jklm /cyc,revers jklm / cyc,revers 

Now the sums are over truly independent moments. Setting coefficients to zero we get the pair 
of conditions R? klm + Ri klm = and R? klm — Ri klm = , whose simultaneous solution is again 
R jkim = o. We conclude that the necessary and sufficient conditions for v 3 to be a symmetry 
of the action are R 3klm = 0. By contracting with the non-singular metric to get a scalar, 

W klm 9j l 9 km = (M-m. (81) 

Since A ^ 2/3, if v is action preserving 

( R jklm = ) } then tr v = and v is automatically 
measure preserving. Thus, non-anomalous linear symmetries of the Yang-Mills model in the 
large- N limit are characterized by those v for which the tensor R) klm vanishes. It suffices to 
check this condition for each word jklm up to cyclic permutations and order reversals. Since 
Rjklm dgpgndg only on the symmetric projection of v 3k , the anti-symmetric part of v jk is 
unconstrained! Thus, a sufficient condition for v %3 to be a non-anomalous symmetry is that it 
be anti-symmetric. However, this is not a necessary condition; there are tracelesq^j v %3 with 
non-trivial symmetric projections for which R) klm = 0. 

Example: We will demonstrate this using the simplest non-trivial example, the 2-matrix 
Yang- Mills model with flat metric gij = 5ij . In this case, the action reads tr S(A) = 



24 tr v = Vj — v jl gij . 
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— 2^ tr L4i, .A2] 2 . Then Vj = v ik Skj- The antisymmetric part of v automatically satisfies 
j^jkim — so us suppose that v l i is a traceless symmetric tensor, i.e. = v 1 ^ and 

v 11 + v 22 = 0. Then the six independent components of R? klm are all identically zero 

i? 1111 = 2{v 11 + v 11 ) - Av 11 = 0, R 2222 = 2{v 22 + v 22 ) - 4v 22 = 0, 

R U12 = 2 ~12 _ ( ~12 + ~12 } = Qj ^1122 = _ ( ~11 + ~2 2) = Q) 

i?i2i2 = 2 {v n + v 22 ) = 0, R 1222 = 2v 12 - (v 12 + v 12 ) = 0. (82) 

So every symmetric traceless v^ k satisfies Ri klm = . We conclude that for A = 2 and g^j = 5ij , 
the Lie algebra of non-anomalous symmetries is s/2 (R) • 

4.3.4 Linear symmetries of Gaussian + Yang-Mills 

For A > 2 let us consider a Gaussian + Yang-Mills matrix model with action 

tr S(A) = X -C^ tr AtAj - ^ tr [A u A,][A k , A x \g ik tf l . (83) 

The simplest case which we will focus on is the two matrix model with flat metric gij = 5ij and 
with covariance a multiple of the identity C l i = m 2 5 l i . In this case the action reads 

tr S(A) = tr ^(Al + A 2 ) - ±- tr [A U A 2 ] 2 . (84) 

We know (sections 14.3. H I4.3.3P that linear non-anomalous symmetries of the Gaussian and 
Yang- Mills parts constitute the o(2) and ^(R) Lie algebras respectively. Their intersection 
is o(2), which is automatically a non-anomalous symmetry algebra of (|84j) . But these must 
be all the linear symmetries, since there can be no cancelation between L v S gauss {G) which 
involves two point correlations and L v Sym(G) which involves 4-point correlations exclusively. 
The corresponding conclusion for A matrix models (again with C %3 a multiple of identity and 
9ij = $ij ) is that the non-anomalous linear symmetries form the orthogonal Lie algebra o(A) . 

4.4 Examples of non-linear non-anomalous symmetries 

We exhibit homogeneous quadratic infinitesimal changes of variable 5Ai = v{ k AjA). which leave 
both action and measure invariant in the large- N limit. In particular, we consider the 2-matrix 
Gaussian with unit covariance, the 3 -matrix Chern Simons model, the 2-matrix commutator- 
squared Yang-Mills model and the 2-matrix Gaussian+YM model. We find a 2,18,6 and 2 
dimensional family of quadratic non-anomalous symmetries in these cases. Moreover, we show 
that quadratic symmetries do not form a Lie algebra by themselves. We demonstrate how to 
obtain non-trivial non-anomalous cubic symmetries via their Lie brackets. 

4.4.1 Quadratic symmetries of the Gaussian model 

Under an infinitesimal quadratic change of variable 5Ai = vj k Ajk , the change in action of a 
Gaussian model with unit covariance is 

5S = tr fiSAiAj = tr S^v^A^ = tr v™ n A mni . (85) 
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Specializing to a 2-matrix model in the large- N limit and taking expectation values, 

L V S{G) = vfGm + {v\ l + v\ 2 + vf)G 112 + {v\ 2 + v\ 2 + vf)G 221 + vfG 222 . (86) 

where we have collected the coefficients of the four independent third rank moments, which are 
all real after accounting for cyclicity and hermiticity. Thus L V S = implies 

11 22 n 11 12 21 22 12 21 f Q n\ 

v{ =v 2 = 0, v 2 = -v[ -Vi, Vy = -v 2 - v 2 . (87) 

To be non-anomalous, v must be volume preserving as well: + vf 1 = 0, which implies (i6~Tj) . 
The solution of this system of linear equations is a two parameter family 

11 22 11 22 n 12 21 12 21 . / 00 \ 

v 1 =v 2 =v 2 = v x = 0, v x = -v x =a, v 2 = -v 2 = b. (88) 

Thus, the non-anomalous quadratic symmetries of the Gaussian model with unit covariance are 

8A 1 = a[A h A 2 ], 5A 2 = b[A 1 ,A 2 ], a,5eR. (89) 

They correspond to the vector fields L Ua b = a(L\ 2 — L 21 ) + b(L 2 2 — L 21 ). The Lie bracket of 
L u b and L Uc d is not a quadratic vector field. Rather, it is a cubic non-anomalous vector field 

[^"Ua,b ' ^«c, J = ( a< ^ ~~ 1 1-^12 ) ^12] ~~ [-^12 ) ^Il] ~~ [-^21 ) ^12] + 1-^21 > ^ll] I 

= {ad-bc){L\ 22 -L\ X2 -2Ll X2 + 2L\ 2X + Ll 21 -Ll u ). (90) 

It corresponds to the one parameter family of infinitesimal changes of variable 

5Ai = (ad-bc)[A 1 ,[A 2 ,A 1 ]}, 5A 2 = (ad - bc)[[A u A 2 ], A 2 ], ad-bc£~R. (91) 

There could, of course, be more cubic non-anomalous symmetries that do not arise as Lie brackets 
of quadratic symmetries. It is satisfying that our point of view tells us something interesting 
even about the Gaussian matrix model. 

4.4.2 Quadratic symmetries of Chern-Simons 

For a homogeneous quadratic change of variable, the change in the expectation value of the 
Chern-Simons action tr S(A) = 2 ^e^ k Aij k in the large- TV" limit is 

L V S(G) = 2iKv l re jki G lm]k . (92) 
To account for cyclicity of Gi m j k we cyclically symmetrize the coefficient. Let 

rplmjk _ y lm e jki _j_ yk^mji , ^jk^lmi _|_ v m 3 £ kli 

Then, L V S(G) = (in/2) £ T lm ^ k G lmjk (93) 

lmjk/~ 

where the sum is restricted equivalence classes of Imjk under cyclic permutations. Accounting 
for hermiticity, we get that L V S(G) = iff 



J2 (T lmjk + T lmjk )$tGi mjk = 0; and £ (T lmjk - T lmjk )%Gi mjk = (94) 

Imjk/r^' Imjk/^' 
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where now the sums are further restricted modulo order reversal. Now we may set the coefficients 
to zero and after adding and subtracting we find 

L V S{G) = o T lmjk = (95) 

where the condition is imposed for all words Imjk modulo cyclic permutations. There are 
c(n = 4, A = 3) = 24 such words for a 3-matrix model. For 9 of these words (1111), (2222), 
(3333), (1112), (1222), (2333), (2223), (3111) and (3331), T lm i k identically vanishes. The 
equations 

T lmjk = for each of the 

remaining 24 — 9 = 15 words are listed below. The words 
are indicated in parenthesis to the left of the equations 

(1212) vf = vf (1123) vf + vf + vf = (3312) vf + vf + vf = 

(1122) vf = vf (1132) vl 1 + vf + vf = (3321) vf + vf + vf = 

(2323) vf = vf (1213) vf + vf - vf - vf = (3132) vf + vf - vf - vf = 

(2233) vf = vf (2231) vf + vf + vf = 

(3131) vf = vf (2213) vf + vf + vf = 

(3311) vf = vf (2321) vf + vf - vf - vf = 0. (96) 

Equations for volume preserving vector fields are vf + v\ % = , or explicitly 

vf + vf + vf + vf + vf + vf = 0, vf + vf + vf + vf + vf + vf = 

and vf + vf + vf + vf + vf + vf = 0. (97) 

We see that (1123) + (1132), (2231) + (2213) and (3321) + (3312) are equivalent to these 
equations. So quadratic symmetries of the CS action are also volume preserving. Moreover 

(1132) - (1123) = (1213); (2213) - (2231) = (2321); (3321) - (3312) = (3132). (98) 

So equations (1213), (2321), (3132) may be discarded. And equations (1122), (2233), (3311) 
are redundant. Thus we are left with 15 — 6 = 9 independent equations for 3 3 = 27 unknowns 
v 1 ™ , leaving an 18 parameter family of non-anomalous quadratic symmetries of the CS model. 



4.4.3 Quadratic symmetries of Yang-Mills model 

The action of the two matrix Yang-Mills model may be written S = — 2~tr[Ai, A2] 2 = ^ tr (vlii22 — 
^1212)- Under an infinitesimal homogeneous quadratic change of variables 6A{ = v{ k Ajk the 
change in the action — ^tr {([5A\, A2] + [Ax, 5A2])[A\, A2]} becomes, in the large- N limit 

- aL v S(G) = vl k L) k (G 12 i2 - G 1122 ) = v?{5^G hjk i, - S{^G hjkl2 ) 
= G U2 i2(2vf + vf + vf) - G U i22(2vf + vf + vf) 

+G 12 i 22 (2vf + vf + vf ) - G ll2 22 (2vf + vf +vf). (99) 

Since the moments that appear are independent, we set the coefficients to zero: 

2vf + vf + vf = and 2vf + vf + vf = 0. (100) 

But these conditions are identical to those for a quadratic vector field to preserve the measure 
of a 2-matrix model (|6ip . Since the above two equations are independent, we have a 2 3 — 2 = 6 
parameter family of non- anomalous homogeneous quadratic symmetries of the 2-matrix Yang- 
Mills model given in (|62p . It is remarkable that every measure preserving linear and quadratic 
vector field also preserves the action of the YM 2-matrix model and CS 3-matrix model in the 
large- N limit. We wonder if this continues to hold for higher rank symmetries or more matrices. 
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Cubic Symmetry: The Lie bracket of two non-anomalous rank- 2 vector fields (if 7^ 0) is a 
rank- 3 non-anomalous vector field (since they form a Lie algebra). This is a way of generating 
new symmetries. Consider two quadratic symmetries of the YM model 

L v = aL 21 ~ 2^11 anc ^ ^ u = ^22- (101) 

which correspond to the choices v^ 1 = a, v\ l = — § while all other vj k vanish, and u 22 = b and 
all other uj k vanish. Their commutator is 

[L v , L u ] = ab ^[L%y, Ll 2 ] — -[L^, Ll 2 ]^ = oh \l\ X2 — h\ 22 + 2^221 + 2^122 j • (102) 

The non-vanishing components of the resulting non-anomalous cubic symmetry are 

w f 2 = ab, wf 1 = ^ab, w\ 22 = ^ab, and wf 2 = -ab. (103) 

One can also check explicitly that this defines a simultaneous symmetry of the action and the 
measure. For example, the change in action is 

StvS = abtv H[A 2 i2 + ~^22i + ~A 122 ,A 2 ] - [Ai, A 222 ]j [A U A 2 ]} = 0. (104) 

The conditions for a homogeneous cubic vector field to be volume preserving are 

wf n + wf m = Q and wf 7171 + wf 1 ™ + w? ni + w? mi = 0, (105) 
and wj kl satisfy these conditions as well. 

4.4.4 Quadratic symmetries of 2-matrix Gaussian + YM model 

Having determined homogeneous quadratic symmetries of Gaussian and YM models, we get 
those for Gaussian+YM model tr S(A) = tr [^(Aj + A$) - ±[A 1 ,A 2 } 2 ] by taking their 
intersection. For, there can be no cancelation between rank 3 & 5 tensors from the action of 
a homogeneous quadratic vector field on the Gaussian and YM terms. Since every quadratic 
measure preserving vector field also preserves the YM action, the intersection is the same family 
SAi = a[Ai,A 2 ], 5A 2 = b[A\, A 2 ],a,b € R as for the Gaussian (section [4.4. lj) . 

5 Supplementing loop equations with Ward identities 

5.1 Gaussian 

LE of the Gaussian are not underdetermined. We do not need the WI in this case. Nevertheless, 
the Gaussian does have non-anomalous symmetries (sections 14.3.11 I4.4.1j) . which lead to non- 
trivial WI. The LE along with WI are an overdetermined system in this case. Nevertheless, 
the WI are consistent with the LE and there is no contradiction, as we have shown in section 
[2j To illustrate this, consider a A -matrix model with unit covariance tr S(A) = \^AiAj . 
The unique solution to the LE states that the odd rank correlations vanish and the even ones 
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Gk are determined by a planar version of Wick's theorem involving a sum over non-crossing 
partitions of K into pairs of indices. For example Gij = Sij, Gijki = SijSki + &iiSjk etc The 
WI corresponding to linear non-anomalous symmetries are (vj are anti-symmetric) 

n n 

L v Gk = v)l{G k = Y, v j 6 i G ki-k^ 1 kk l+1 ...k n = J2 v ki G ki-h- 1 kk l+1 ...k n =0. (106) 

i=l i=l 

Thus, the moments must be o(A) invariant tensors. We can check that these WI are consistent 
with the LE. For example, with n = 1 we get the WI v\G k = 0, for all antisymmetric vf . But 
there are anti-symmetric vi with non- vanishing determinant, and Gj must lie in their kernel 
which is trivial. So Gj = 0, as implied by the LE. The WI for odd n are trivially satisfied by 
solutions to the LE, since odd rank moments vanish. WI for n = 2 are 

vt 1 G k k2+vt 2 G kl k = Q, V antisymmetric vf. (107) 

For Gij = , the LHS becomes + which vanishes on account of anti-symmetry of v . 
Similarly we can check that the WI are consistent with the LE for n = 4, 6, • • • . 

We could do the same for quadratic symmetries. Let us consider the two-matrix Gaussian 
model. The WI L v Gk = following from quadratic non-anomalous symmetries (section 14.4. ip 
L v = a{L\ 2 — L\i) + b{L\ 2 — for arbitrary real a, b are 

Sk 1M G L i2M = and 5f M G L21M = 0. (108) 

These WI are consistent with the LE (for \K\ < 4, that we checked, these WI are consequences 
of cyclicity and do not contain new information). For more nontrivial use of the WI we must 
progress to non-Gaussian multi-matrix models whose LE are underdetermined. 



5.2 Gaussian plus Yang-Mills 

The matrix integrals for correlations of the 2 -matrix Gaussian + YM model, whose action is 

tvS(A)= tT[^(Al + Al)-±[A u A 2 n (109) 

converge. Recall that the commutator of hermitian matrices is anti-hermitian, and the square of 
an anti-hermitian matrix is non-positive. Thus, the quartic term is non-negative. The quadratic 
term ensures that as any matrix element goes to ±oo, the action goes to +oo. Thus, the 
Boltzmann weight e~ N tr s vanishes at least exponentially fast as any matrix element goes to 
±oo. Thus, all polynomial observables have finite expectation values. From this we conclude 
that the LE and WI are rigorously valid. In section [231 we obtained the LE for |/| < 4. They 
left a number of correlations undetermined. In section [4.3.4l we found that linear non-anomalous 
symmetries of this model form the o(2) Lie algebra parameterized by v *■ such that v\ = v 2 = 
and v 2 = —vf . The corresponding WI, which we will use to supplement the LE, read TGk = 
for all words K , where T = L\ — L 2 . These are listed in appendix [B] for moments of rank up to 
4. They imply that all Gi vanish. The only Gij that might be non- vanishing are G\\ = G22- 
All 3-point Gijk vanish. 4-point correlations vanish except possibly Cam, G2222, G1212, G1122 
and their cyclic permutations. They must, however, satisfy the relations G\\\\ = G2222 and 
Gini = 2G1122 + Gi2i2- Some of these conditions could also have been got from the LE, (|2.2p . 
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We need one more condition on rank- 2 moments and two more conditions on rank- 4 moments 
to determine all moments of rank < 4. The LE for |/| = 1 gives one new condition 

Gii = l + -(2G 12 i2-2Gii 22 ). (110) 
a 

The LE for |/| = 2 fsection l2.2p relate 3 and 5 point correlations. Using the fact that all 3-point 
correlations vanish, they tell us that G11212 = G11122 and G12122 = G11222 • Supplementing these 
LE with the WI for 5-point correlations TGij^im = V ijklm (which we do not list explicitly), 
we are able to conclude that all rank-5 correlations vanish. 

Thus far, we have found that the only correlations with rank < 5 that could be non- 
vanishing are Gu, G22, Gnu, G2222, G1212 and G1122 > up to cyclic permutations. We have found 
4 relations among these 6 unknowns: 

G11 = G22, Gnu = G2222, Gini = 2G1122 + G1212, G11 = H — (2G1212 — 2G1122) (HI) 

a 

Thus, by use of the WI, we have reduced the underdeterminacy of the LE. We could proceed 
further in this manner. The LE for |/| = 3 relate rank-4 and rank-6 moments, while the WI for 
rank- 6 moments give further conditions on rank-6 moments. We could also look for additional 
conditions using the WI from quadratic symmetries found in section 14.4.41 Du t we postpone that. 
Our purpose here was only to illustrate the general framework we have developed. In a separate 
paper, we hope to return to a more thorough study of the correlations of this model using the 
LE and WI and their comparison with other approaches [26\ I27j or monte-carlo simulations |28|. 

5.3 Chern-Simons 3-matrix model 

The CS 3-matrix model has action tr S(A) = ^e ijk tr A ijk = 2in tr AiL4 2 ,A 3 ]. We expect 
its matrix integrals to diverge. To see this, go to a basis where A<± is diagonal, then the action 
is independent of the diagonal elements of ^3 , due to the commutator. So integration over 
the diagonal elements of ^3 would diverge. Our derivation of the LE and WI holds at best 
formally for this model. We do not know whether the LE and WI are a consistent system for 
this action. Nevertheless, we consider them formally to illustrate the general framework. We 
find no inconsistency, at least for correlations up to rank 3 . The LE of the CS 3-matrix model 
are under determined (section I2.2j> . WI corresponding to non-anomalous linear symmetries (I75p 
are obtained by imposing the conditions L v Gk = for all K and traceless Vj : 

v{S§ liK2 G KljK2 = V traceless v). (112) 

These are the conditions that Gk be (cyclic and hermitian) invariant tensors of SLs(R). We 
will work out the WI explicitly for \K\ = 0,1,2,3. For K empty, this is a vacuous condition, 
so put K = k to get the WI viGj = for all traceless vi. There are traceless v k with non- 
vanishing determinant, and Gj must lie in the kernel of such linear transformations. But this 
kernel is trivial, so Gj = 0. For K = kl we get 

vjGkj + v J k Gji = V traceless v. (113) 

Putting k = I = 1 we get v\G\\ + i>iGi2 + V1G13 = 0. But this must hold for all real v\,v{,v\ 
so that G11 = G12 = G13 = 0. Putting k = I = 2 we get G21 = G22 = G23 = 0. Finally putting 



28 



k = I = 3 we get v\G^\ + — (v\ + Again U3, vf, i;*, uf are freely specifiable so that 

G31 = G32 = G33 = 0. From this we conclude that all Gy = 0. This is of course consistent 
with the remaining WI gotten by putting k = 1, 1 = 2 etc since Gij = is an obvious solution 
of the homogeneous system vjGkj + v^Gji = 0. 

As for wffil for rank-3 correlations^], we set K = Mm and get 

v k G jim + viGkjm + v J m G k ij = V traceless v. (115) 

k = I = m = 1 gives v\G\u + v\Gi\2 +ufGii3 = whence Gm = G112 = G113 = 0. Similarly, 
putting k = I = m = 2 and = / = m = 3 we get G122 = G222 = G223 = G133 = G233 = G333 = 
. The only remaining undetermined rank-3 correlations are G123 and G132 . The remaining 
WI are either vacuous (e.g. k = 1, 1 = 2, m = 3) on account of v being traceless or (e.g. 
k = l = i )m = 2) give G123 + G132 = 0. 

To summarize, the WI due to linear non-anomalous symmetries imply that all correlations 
of rank < 3 vanish except for G123 and G132 (and their cyclic permutations), and these are 
related by the WI G123 + G132 = . WI remedy the underdeterminacy of LE of the CS model. 
The only non-trivial condition from the LE was (see section 2.2.2 of ref . |20| ) 

9(G 12 3 " G132) = (116) 

This, along with the WI G123 + G132 = now allows us to determine all correlations up to rank 
3, the only non- vanishing ones (up to cyclic symmetry) ara^l 

G123 = -77- and G132 = -— (117) 

We checked that this result is consistent with the WI corresponding to quadratic non-anomalous 
symmetries obtained in section 14.4.21 So at least up to rank-3 moments, the WI cure the under- 
determinacy problem of the LE! We could proceed in this manner to higher rank correlations. 



5.4 2-matrix Yang-Mills: A cautionary tale 

The matrix integrals for correlations of the YM 2-matrix model tr S(A) = — 5— tr [Ax, A2} 2 do 
not converge }25j due to a similar argument as given for the CS (section l5.3p . Thus, our derivation 
of the WI and LE is not strictly valid. We cannot be certain that they form a consistent system. 
In fact, we find that the WI and LE for this model do not form a consistent system when 
considering rank-4 correlations. Despite several checks, we could find no calculational error. We 
do not know the deeper reason for this inconsistency, but suspect it could have something to do 
with the lack of convergence of matrix integrals invalidating our derivation of the WI and LE. 
Thus, it is probably good to be cautious in formal use of the WI and LE. 

25 LI (no sum over i) is a number operator, it counts the number of i's in a given moment. For example, 
consider GujuciL , where none of the multi-indices /, J, K, L contain an i . Then 

L\GliJiKiL ~ GliJiKiL + GliJiKiL + GujiKiL = SGnjiKiL- (H4) 

The number operators commute [Z4,L^-] = 0. Thus V_\ L\ measures the rank of a given moment. This can be 
used to get most of the rank-3 moments of the CS model by employing WI involving L\ — L% and L\ — Li, . 

26 Accounting for cyclicity and hermiticity, the space of rank-3 tensors is 11 dimensional, see appendix iDl 
27 Hermiticity and cyclicity mean G123 — G132 which implies KG123 = KG132 and ^sGwa = ~ 9Gi32 
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The LE of the YM 2-matrix model are under deter mined, (section I2.2p . Recall that the LE 
do not determine any moments of rank 1,2 or 3. Here, the WI come to the rescue. Recall 
(|4.3.3|) . that the non-anomalous linear symmetries of this model form the Lie algebra s^R), 
spanned by L\ — L 2 ,, L\, L\ ■ The WI L v Gk = for \K\ < 3 suffice to determine all 1, 2, and 3 
point correlations, and imply they are all zero. Let us also consider WI for moments of rank 4: 

{L\ — -^DGiiii = 4Giiii = 0, (L\ — LDG2222 = — 4G2222 = 0, 
(L 1 — L 2 )Gni2 = 2Gni2 = 0, L 2 Gni2 = G1212 + 2G1122 = 0, 
-^2^1212 = 2G2221 = => G2221 = 0. (118) 

While the WI determine many rank-4 correlations, they give only one relation between G1212 
and G\\22 ■ To fix them we use the simplest conditions coming from the first of the two LE 

1=1: 2G1212 — 2G1221 = —a ; 1 = 2: 2G2212 — 2G2221 = 0- (H9) 

The second equation is vacuous, but thanks to the WI we see G1222 = 0. As for the first 
equation, WI provides us with the condition G1212 + 2G1122 = . So we get 

a, a 
Gm2 = - T and G 1 i22 = 7r - (120) 
3 o 

However, this is not consistent with the WI for non-anomalous quadratic symmetries (section 
14.4. 3p . In particular the WI obtained from linear symmetries together with 

1 1 

{L 12 — 2^22)^112 = G1212 + -G1122 = (121) 

implies G1212 = G1122 = , which is a contradiction. 



6 Some outstanding questions 



A summary and discussion of the results of this paper was given in the introduction. Here, we 
list some questions raised by our work. (1) We have only addressed the exact determination 
of normalized correlations in large- N matrix models using the LE and WI. But what about 
the partition function or free energy? (2) It is interesting to know whether the LE and WI 
together determine all single-trace correlations in the large- N limit. (3) We have only discussed 
infinitesimal non- anomalous symmetries. Many models also possess discrete non- anomalous 
symmetries, which lead to useful relations among correlations. Some of these relations are 
actually a consequence of the LE or WI. But in general, it may be necessary to supplement 
the LE and WI by conditions from discrete symmetries. (4) Detailed study of LE and WI of 
specific multi-matrix models should clarify whether we need additional conditions to solve for 
all correlations. (5) It is interesting to identify matrix models with a maximal family of non- 
anomalous symmetries. Interestingly, we found that the 3-matrix CS model and the 2-matrix 
commutator-squared YM models each possesses a maximal family of linear and quadratic non- 
anomalous symmetries. (6) It is interesting to classify the solutions to the simplest WI. For 
example, the correlations that satisfy WI for linear symmetries of the Gauss+YM model must 
be invariant cyclic hermitian tensors of the orthogonal Lie algebra. What is the general form 
of such tensors? (7) We observed that for n > 1 , Lie brackets of rank n non-anomalous vector 
fields are rank n + 1 non-anomalous vector fields, provided they are non- vanishing. It would 
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be interesting to study this Lie algebra of non-anomalous symmetries in specific examples. Can 
it be infinite dimensional? If so, might the model be integrable in some sense? (8) We wonder 
whether the full gauge fixed Yang- Mills theory in the large- N limit has any additional non- 
anomalous symmetries besides Poincare invariance and BRST invariance. Our work indicates 
that such symmetries can be far from obvious and highly non-linear. 

Acknowledgements 

We thank Jean Yves Thibon for discussions on the dimension of the space of cyclically symmetric 
tensors. We also thank S. G. Rajeev, G. Arutyunov, S. Vandoren O. T. Turgut and G. 't Hooft 
for discussions. GSK thanks the Feza Gursey Institute, Istanbul for hospitality while a part of 
this work was done, and support of the European Union in the form of a Marie Curie Fellowship. 

A Alternative derivation of SDE preserving hermiticity 

Consider the change of variables corresponding to translation by a constant Hermitian matrix e 
(A t ) a b -> {Ai)l + el, and - (A-fi for j^i (122) 

in the integral ( [dA] = J[ k dA^ is the Lebesgue measure on independent matrix elements) 

/= f[dA] (A/)g <$> K e - NtrSJAj . (123) 

Here <&k = jjAk ■ The value of the integral should be unaltered under this change of variables^] 

SI = e d J [dA] (A h ) a c (A l2 ) d b <J? X e 

+ el J [dA] [AtK S% m ±(A ML ) d c e 

- e c d J [dA] (Ajy b $ K NS J 5 J / J * (Aj 2 . h ) d c e ~ NtlSJ ^ = (124) 

where we also used the translation invariance of the measure. Since this holds for arbitrary 
Hermitian e we conclude § 

S^HAjXiAi^K) + p L K tM ((Aj) a b (A ML ) d c ) = NS^dA^KiAj^i). (125) 
Contracting a with c and b with d, and dividing both sides by N 2 we get 

S^^^j^k) + 6% M ±($ LIM ) = S w *($ JlIJ2 $ K ). (126) 
Since this must hold for every / and i , it is equivalent to the equations 

vIS^^^^^k) + vl5^ M ^um) = vIS^^^j.jj^k) V v\ G R. (127) 

In a completely analogous fashion, we repeat this calculation with several insertions • • • &k u 
and get the Schwinger-Dyson equations obtained earlier in ([7|). 



-N tr S J Aj 



-N tr S J Aj 



2S The following formula is useful: = S 1 /" 2 (AiJlSfS^AiJl = S' 1 " 2 (A h )%{Ai 3 )* . Contracting a with 

b one gets |^ = 6?"' (A^A^t = 6?"* {A Wl ) d c . 

29 It is shown in Adler's book [29] page. 26 that if tr (eH) —0 for all Hermitian e then H = 
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A.l Other possible changes of variables in matrix integrals 



LE are under determined, so do the Gi satisfy other equations? The WI are such equations and 
with the LE, may go a long way towards fixing Gj . Here we consider two other types of changes 
of variable in matrix integrals to see if they give new equations. However, we do not find any. 

Consider an infinitesimal change SAi = v(Aj where v\ is a Hermitian matrix for each / & 
i; previously they were real numbers. Conditions for invariance of partition function are 

- NS hlh tr (A Jhh vi) + 5 J / jJ * tr (vjAjJ tr A, h = (128) 

Since these must hold for arbitrary matrix elements [v{]% we get 

- NS hU *(A Jhh )t + Sj^iAjXb tr A, h = (129) 

But these are not equations for trace invariants. To get equations for Gj , we must take a trace, 
divide by N 2 and take expectation values. But this leads to the LE derived before. 

Next we consider an infinitesimal change SAi = [v(,Ai] with arbitrary hermitian matrices 
v{ . These types of change of variable do not appear in BRST transformations but do appear in 
gauge transformations A^ — > A^ + d^A + [A^x), A(x)] . Invariance of partition function implies 

5j ljh tr (vjA h ) tr A h - 5j ljh tr A h tr (vj A h ) + S hjh [A hh , A/] = 0. (130) 

Using the arbitrariness of v{ we get 

5^ h (A h r b tr A h - S T / jh tr A h (Ai 2 )i = S h ^[A hh , Aj} a b . (131) 

As before, we must take a trace to get an equation for the G\ , but in fact we get a triviality. 
Thus, we have not found any equations for Gj in addition to the LE and WI. 



B WI for Gaussian+YM model 



Below is a list of Ward identities TGr = for moments of rank up to \K\ = 4 in the 2 -matrix 
Gaussian+YM model. They correspond to the non-anomalous linear vector field T = L\ — L\ . 

TGi = -G 2 = 0; TG 2 = G 1 = 0; TG U = -G 2 i - G 12 = -2G 12 = => G 12 = 0; 

TG\2 = Gn — G22 = ^> G11 = G22; TG 22 = 2G12 = G12 = 0; 

TGm = — 3Gn2 = => Gu 2 = 0; TG\\ 2 = Gm — 2G122 = => Gm = 2G122 = 0; 

TG122 = 2G112 — G222 = G222 = 2G211 = 0; TG222 = 3Gi22 =y~ G\22 = 0; 

TG\\\\ = — 4Gni2 = 0; TGim = Gim — 2G1122 — G1212 = Gim = 2G1122 + G1212; 

TG1122 = — 2G1222 + 2Gni2 = G1112 = G1222; TG1212 = 2Gni2 — 2G2221 = 0; 

TG1222 = 2G1122 + G1212 — G2222 = 0; and TG2222 = 4Gi222 = 0. (132) 



C Is there a model for which Ward identities determine all G{1 

Are there any non-trivial multi-matrix models where the WI are sufficiently numerous to de- 
termine all (or a maximal set of) correlations without using the LE? We seek models with a 
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very large family of non-anomalous symmetries. Trying to answer the corresponding question in 
two dimensional quantum field theory has proven very fruitful, as evidenced by the progress in 
2d conformal field theory. The latter are so symmetrical that a maximal family of correlations 
can be determined by milking conformal invariance. Here we make an elementary observation. 
The WI L v Gk = are a system of homogeneous linear equations. So they are either un- 
derdetermined (if the determinant of the system vanishes) or admit only the trivial solution 
Gk = 0, MK . Though the WI can give us much information on correlations, they cannot deter- 
mine all of them except in the trivial case where they are all zero. For example in a 2-matrix 
model, if we consider an extreme (and probably unrealistic) case where all measure preserving 
vector fields are also action preserving 30 !, then it follows that all correlations of rank up to 4 
vanish. This leaves open the question of identifying non-trivial models with a maximal family 
of non-anomalous symmetries, i.e. the ones for which the WI are most useful. 



D Cyclically symmetric tensors of rank n 



The real dimension of the space of cyclic hermitian tensors G f i 1 ...j n on a vector space V of 
dimension A (i.e. 1 < i±, ■ ■ ■ , i n < A) is 

c(n,A) = -J20(d)A- n/d (133) 

a\n 

where 4>(d) is Euler's totient phi-functional and the sum is over all divisors of n. c(n, A) is the 
number of independent correlations of rank-n in a A matrix model. We thank the mathematician 
Jean Yves Thibon of Universite de Marne-la-Vallee, France for sharing this formula with us. This 
answers a question posed in Appendix A of Ref.[20j. It comes from the character of GL(V) 

1 A 

— (p(d) p^ d (S,) where Pd(C) = ^2^i * s P ower sum symmetric function. (134) 

n d\n i=l 
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